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Abstract 

A nested Fermi surface with nearly parallel orbit segments is found to yield 
a singlet d-wave superconducting state at high temperatures for a restricted 
range of the on-site Coulomb repulsion that avoids the competing spin den- 
sity wave instability. The computed superconducting transition temperature 
drops dramatically as the nesting vector is decreased, in accord with recent 
photo-emission data on the Bi2212 and Bi2201 cuprates. 
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Even though the BCS theory jl] provides a successful description of conventional su- 
perconductors, the concept of alternate electron pairing states of finite angular momentum 
has evolved to encompass various physical systems. Spin fluctuations suppress 0,|J the 
BCS isotropic pair binding induced by phonon exchange, and hence materials with strong 
repulsive interactions are prospects for anisotropic pairing. Superfluid He 3 exhibits p-wave 
pairing [||], and heavy Fermion superconductors offer another unconventional case at very 
low temperatures. 

Copper oxides with high superconducting transition temperatures T c exhibit abnormal 
electrical transport and optical properties. The observed linear temperature variation of the 
resistivity was attributed by Lee and Read to electron-electron collisions on a perfectly 
nested Fermi surface in the form of a square, and the anomalous linear frequency variation 
of the damping has been derived for a partly nested Fermi liquid (NFL). |7]|| However, 
nesting models need to consider the competing spin density wave (SDW) instability. 

The evident correlation of high T c values in cuprates with departures from the standard 
Fermi liquid behavior suggests that the physical origin of the anomalous damping is a key 
source of the superconductivity. This connection is found here by the present microscopic 
theory for a partially nested Fermi surface. The on-site Coulomb repulsion U provides the 
primary interaction, while the nested orbit topology is the key determinant of the attraction. 
The pairing symmetry and binding via exchange of spin fluctuations are determined by the 
nesting vector Q. The d angular momentum state is found to be favored in the Bi2212 
cuprate where our calculations reveal the existence of superconductivity for values of U that 
avoid the SDW instability. 

Our investigation was inspired by the discovery of Scalapino et al. || that d-wave pairing 
may be possible in a Hubbard model. Similar correlations were found in Monte Carlo 
simulations [[TIJ on a small lattice. However, for the band fillings and Fermi surface topologies 
treated by Scalapino and others [PH-|T5fl, the lowest order estimates give very small T c 
values. Other anisotropic pairing proposals have been applied to organic compounds |TJ 
and heavy Fermion superconductors fl5| . 



We consider the Hubbard Hamiltonian 

where the electron (or hole) energy band E(k) is represented by the tight-binding expression 

E(k) = — 2t(cos k x + cos k y — B cos k x cos fc y + — ), (2) 

Z7 is the Coulomb repulsion between electrons at a given site, and ct (cg CT ) represent creation 
(destruction) operators of momentum fc and spin cr. The Fermi surfaces for this model are 
shown in Figure [l] for a rounded orbit (Fermi liquid FL), and a nested surface that resembles 



photoemission experimental data by Dessau et al. [16j and Shen et al. [[17]] on E^S^CaC^Og. 

Electron scattering in the singlet spin channel involves a direct Coulomb term and the 
exchange of spin fluctuations |4j] of the form shown in Figure 2. The phase space for the 
scattering and the susceptibility enhancement near the nesting vector are important for the 
d-wave pairing as well as for the SDW instability condition Ux' = 1- 

We compute the real part of the susceptibility, x'(g, tu), using the standard definition 
]Tj| and include self energy corrections of the NFL form J7J. The results for the E(k) model 
chosen to represent Bi2212 are shown as a function of momentum in Figure 3. The double 
peak structure at low frequencies is similar to the neutron spectra for the La2- a; Sr :r Cu04 
superconductor |19| . Another consequence of nesting is the scaling of the spin susceptibility 
as a function of uj/T which has been confirmed by neutron scattering on several cuprates 
2"0fl . The SDW constraint on the susceptibility requires U < 1.1 eV in the case of Bi2212 
where the bandwidth is estimated to be 1.5 eV from photoemission data. 

Decomposition of the two-particle scattering into angular momentum channels yields an 

effective pairing coupling 

A = Sjgfr 9iik) V(k,k') gjp) S[E(k)] 5[E(k')} 
Z- k 9f(k) S[E(k)] 



The conventional symmetry classification of the basis set g\ f9|-|T3"|| is g s — 1 for the s-wave 
states, g p = sink x for p- waves and the d-wave states with g x 2_ y 2 = cosk x — cosk y and 
9xy = si n k x sin k y . 



The primary pairing interaction V(k, k') (see Fig. is given by the term with two spin 
fluctuation bubbles, U 3 x' 2 (k ~ k'), and the exchange term proportional to U 2 x'{k + k'). If 
X'(q, 0) is approximately constant, an inspection of Eq. [5] reveals that superconductivity by 
this mechanism is not possible because A < 0. 

To calculate the coupling, we represent the susceptibility by a Gaussian form 



X \q,0) = A + B(Q)ex V 



(\q*\ - Qx) 2 + (\q y \ - Qy? 



(4) 



2a 2 

where q = k — k',Q = (n,7r), and the constants A and B(Q) determine the normalization 
for the Gaussian. This model yields a reasonable fit to the computed Bi2212 susceptibility 
shown in Figure |3|. The actual nesting peaks for our Bi2212 model are at Q\ = (£7r, it), 
Q2 = (k^tt), Q3 = (27T — ^7r,7r) and Q± = (vr,27r — £7r), with £ = 0.91, but this four-peak 
structure produces only small corrections to the T c determined by the simple Gaussian in 
Eq. § 

Our previous analytic derivation J7| of the NFL susceptibility using the nesting approxi- 
mation E(k+Q) = —E(k) gave a logarithmic temperature variation of Xnfl(Q> 0)- However, 
the present calculation gives a smaller susceptibility (see Fig. ||), with a weaker temperature 
dependence. The susceptibility reduction is caused by the rounded corners in our Fermi 
surface model and the influence of the NFL self-energy V = Max(T, \u\). Together, these 
features avoid the SDW formation at intermediate values of the interaction, such as U ~ 0.96 
eV (compared to the bandwidth 8t = 1.5 eV) that is used here. 

Numerical integration over momenta gives the coupling X x 2_ y 2, and our computed sus- 
ceptibility indicates an energy cut-off u c ~ 0.3 eV. The leading order evaluation of the 
superconducting temperature becomes 

T c = u; c exp(- — —J. (5) 

\A x 2_ y 2 ) 

We first find that the Fermi Liquid topology (dashed curve in Figure |1|) gives \ x 2_ y 2 = 
0.016 which corresponds to a vanishing T c . This example is qualitatively similar to other 
cases studied by several groups [PH^3]|. The random phase approximation (RPA) contribu- 



tions enhance this coupling near the SDW instability ||[Ij]-|l3[] . 



Nesting topologies increase the attraction in a d-wave channel as we demonstrate for 
the Bi2212 cuprate model. Using band parameters that yield the nested Fermi surface in 
Fig. [1|, U = 1.0 eV and x'max = 0-92 eV _1 , we obtain a d-wave coupling in the lowest order 
\ x 2- y 2 = 0.27 that gives T c = 90 K. 

If the Coulomb coupling is of intermediate strength, e.g. U\' = 0.9, the RPA enhance- 
ment of the nested case would elevate the coupling to A x 2_ y 2 = 2.7 and thereby predict an 
enormous T c . This situation should stimulate further research on vertex corrections and self 
energy effects that may offset the RPA series enhancement. 

The sensitivity of the coupling to the magnitude of the nesting vector is illustrated in 
Figure (|. The Bi2201 cuprate exhibits a low T c = 6 K despite having abnormal resistivities 



and optical properties in league with the high T c cuprates. Shen et al. |21| have discovered 
by photoemission spectroscopy that this cuprate possesses a nesting vector close to 0.8Q, as 
compared to the ideal half filled case of Q = (jr, 7r) and the Bi2212 situation with a nesting 
vector of 0.9Q. This correlation is compatible with a sharp drop in the calculated T c values 
as seen in Figure [|. Our model predicts that the spin susceptibility peaks seen in Figure |3| 
for the Bi2212 case should spread apart in Bi2201 and this feature may be tested by neutron 
scattering measurements. 

We do not find superconductivity of xy symmetry for our Fermi surface geometry. The 
x 2 — y 2 state for the present nesting model is consistent with photoemission measurements of 



the energy gap anisotropy in Bi2212 [17]. If the Fermi surface is rotated in other cuprates, as 



suggested by photoemission spectra of YBa2Cu3C , 7_ < 5 by Liu et al. [p2j] , then states of other 
symmetry should be examined in more detail. We find a vanishing T c for p-wave symmetry 
pairing using g x = sin k x in both the Fermi liquid model and the nested Fermi surface. 

Our model for Bi2212 yields a Van Hove singularity in the density of states that is located 
0.04 eV~ 500 K above the Fermi energy. An arbitrary increase of the chemical potential 
towards the logarithmic singularity would result in a SDW phase as shown in Fig. ^ 

Impurity scattering should be detrimental to anisotropic pairing as well as for the SDW 
formation. By analogy with the Abrikosov-Gorkov theory, |23[ d-wave suppression by dis- 



order constrains T c in the cuprates |24]]. Similarly, non-magnetic impurities also impede 
the competing SDW transition P5| . The case of chromium reveals a further sensitivity of 
the SDW to impurity induced shifts of the band structure ||26|| . Impurities at sites in the 



copper oxide planes should be more destructive for d-wave superconductivity than those at 
interplanar sites. 



The origin of nesting features in cuprates is evident in band structure calculations |27 
because of the nearly half-filled d-bands in two dimensions. Logically, the relative persistence 
of parallel segments in a given band subjected to doping may be stabilized by a second band 
that acts as a charge reservoir. 

Theoretical extensions of the present work may be relevant to higher order spin fluctua- 
tion graphs, including the "spin bag" variety ||28|| , and the self-energy and vertex corrections. 
Nesting of a two-dimensional electronic structure produces |7||| a linear frequency variation 
of the quasiparticle damping that bears similarities to the Luttinger theory J29| for a one 



dimensional electron gas, which also exhibits remarkable charge and spin dynamics |3(| , pT 

Nevertheless, nesting in two dimensions is distinguished by a crossover temperature T* 
below which the electronic response reverts to standard Fermi liquid behavior. Accordingly, 
the concept of a well-defined Fermi surface is valid in the NFL approach, despite the unusual 
damping features that arise above T* and a corresponding frequency crossover uo* that are 
determined by the nesting geometry. 

Our analysis may provide a guide to the design of new superconducting materials. The 
primary ingredients for the d-wave pairing are a Fermi surface topology with a nesting 
vector restricted to a narrow range, a Coulomb repulsion U of intermediate strength, and a 
planar electronic structure that accepts intercalant atoms with weak impurity scattering of 
electrons (or holes) in the conducting layers. 

We have benefited from discussions with J. P. Collman, D.S. Dessau, D. Huse, R.B. 
Laughlin, W.A. Little, and Z.X. Shen. We (J.R. and C.R.) appreciate the hospitality of the 
physics department at Stanford University during a visit sponsored by a Sesquicentennial 
Associate award from the University of Virginia. 
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FIGURES 

FIG. 1. A nested Fermi surface (NFL) shown by the solid curve was calculated to fit the 
photoemission data points of Dessau et al. [16] using the tight-binding model of Eq. 2 with 
B = 0.165 and \i = —0.56. The nesting vector is Q* ~ 0.91(7r, 7r) in this case. By contrast, the 
dashed curve for the same value of B but a chemical potential ji = —1.6 shows a rounded orbit 
reminiscent of a standard Fermi Liquid (FL). 

FIG. 2. Diagrams for the electron-electron scattering in the singlet channel show the direct 
Coulomb repulsion by a dotted line and the spin fluctuation exchange processes with a bubble 
representing the susceptibility. In the d-wave channel for a nested Fermi surface, the leading order 
attractive contributions from the graphs involving the susceptibility are of the same order, whereas 
the direct bare Coulomb repulsion gives no contribution in the Hubbard model, because U is 
assumed to be momentum-independent. 

FIG. 3. The calculated real part of the susceptibility for the Bi2212 band parameters B = 0.165 
and /U = —0.56 is shown as a function of momentum \q\ along the direction q x = q y by the solid 
curve for the damping Tnfl = Max(T, |u;|), and by the dashed curve for a damping T = 0. The 
calculated maximum x'(Q*) — 0-97 eV _1 constrains U < 1 eV which compares with the bandwidth 
8t = 1.5 eV that we estimated from the photoemission data of Ref. 16. The dot-dashed curve 
represents the Gaussian model of Eq. 4 with a = 1.2. 

FIG. 4. T c values as a function of the nesting vector Q* calculated using the tight-binding model 
of Eq. 2 are shown by the solid curve. The Bi2212 parameters are B = 0.165 and \i = —0.56. The 
points are the experimental values. A surface with a nesting vector Q* ~ 0.8(7r, it) appropriate to 
the Bi2201 photoemission data [22] was simulated using B = 0.33 and fi = —1.36 which lowers T c . 
Intermediate Q* cases were found by linear interpolation of the band structure. The shaded region 
designates the spin-density wave (SDW) regime. 
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